The Strong CP Problem and Axions 



R.D. Peccei 

Department of Physics and Astronomy, UCLA, Los Angeles, California, 90095 
peccei@physics .ucla. edu 



Summary. I describe how the QCD vacuum structure, necessary to resolve the 
U(1)a problem, predicts the presence of a P, T and CP violating term proportional 
to the vacuum angle 9. To agree with experimental bounds, however, this parameter 
must be very small (9 < 10 -9 ). After briefly discussing some possible other solutions 
to this, so-called, strong CP problem, I concentrate on the chiral solution proposed 
by Peccei and Quinn which has associated with it a light pseudoscalar particle, the 
axion. I discuss in detail the properties and dynamics of axions, focusing particularly 
on invisible axion models where axions are very light, very weakly coupled and very 
long-lived. Astrophysical and cosmological bounds on invisible axions are also briefly 
touched upon. 



1 The U(1)a Problem and Its Resolution 

In the 1970's the strong interactions had a puzzling problem, which became 
particularly clear with the development of QCD. The QCD Lagrangian for 
N flavors in the limit of vanishing quark masses to/ — > has a large global 
symmetry: U(N)v x U(N)a- Since m u , rrid << Aqcd, we know that, at least 
for these quarks, the limit of sending the quark masses to zero is sensible. Thus 
one would expect the strong interactions to be approximately U(2)v X U(2)a 
invariant. 

What one finds experimentally is that, indeed, the vector symmetry cor- 
responding to isospin times baryon number [U(2)v = SU(2)i x U(1)b] 
is a good approximate symmetry of nature, as manifested by the appear- 
ance of nucleon and pion multiplets in the spectrum of hadrons. For axial 
symmetries, however, things are different. Dynamically, quark condensates 
< uu >=< dd >^ form, breaking the axial symmetry down spontaneously. 
So one does not expect approximate mixed parity multiplets in the hadronic 
spectrum, but vestiges of the four Nambu-Goldstone bosons associated with 
the breakdown of U(2)a- Although pions arc light, ~ 0, there are no signs 
of another light state in the hadronic spectrum, since mi >> m\. Weinberg 
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dubbed this the U(1)a problem and suggested that, somehow, there was 
no U(1)a symmetry in the strong interactions. 

The resolution of the U(1)a problem came through the realization by 
't Hooft [2] that the QCD vacuum has a more complicated structure. The 
more complex nature of the QCD vacuum, in effect, makes U(1)a not a true 
symmetry of QCD, even though it is an apparent symmetry of the QCD 
Lagrangian in the limit of vanishing quark masses. However, associated with 
this more complicated QCD vacuum there is a phase parameter 6 and only 
if this parameter is very small is CP not very badly broken in the strong 
interactions. So the solution of the U(1)a problem begets a different problem: 
why is CP not badly broken in QCD? This is known as the strong CP problem. 

A possible resolution of the U(1)a problem seems to be provided by the 
chiral anomaly for axial currents. 3 The divergence of the axial current 
J£ associated with U(1)a gets quantum corrections from the triangle graph 
which connects it to two gluon fields with quarks going around the loop. This 
anomaly gives a non-zero divergence for J§ 
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where F aflu = \t^ va ^F^ . Hence, in the massless quark limit, although for- 
mally QCD is invariant under a U(1)a transformation 

qf e »75/ 2(?/ (2) 
the chiral anomaly affects the action: 

SW = aj d A xd^ = J d 4 xFrP a ^. (3) 

However, matters are not that simple! The pseudoscalar density entering 
in the anomaly is, in fact, a total divergence 0|: 

F£»F alu , = d»K», (4) 

where 

K" = e^A aa [F af}l - 9 -f abc A bfi A cl ] . (5) 
Because of these identities SW is a pure surface integral 

SW = J d 4 xd^ = J dv M K". (6) 

Hence, using the naive boundary condition that A£ — at spatial infinity, 
one has / da^K^ = and U(1)a appears to be a symmetry again. What 't 
Hooft 2 showed, however, is that the correct boundary condition to use is 
that A^ should be a pure gauge field at spatial infinity (i.e. either A£ — 0, or a 
gauge transformation of 0). It turns out that, with these boundary conditions 
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there are gauge configurations for which J da^K^ ^ and thus U(\)a is not 
a symmetry of QCD. 

This is most easily understood by working in the A° gauge. Studying 
SU(2) QCD for simplicity, in this gauge [5] one has only spatial gauge fields 
A l a . Under a gauge transformation these fields transform as: 

^ Ta A l a = A 1 ^ QA'Q- 1 + -\7 i QQ- 1 . (7) 

Thus vacuum configurations either vanish or have the form -V 1 ]?]? -1 In the 

9 

A° = gauge one can further classify these vacuum configurations by how J? 
goes to unity as r — * oo: 

[2 n ^e a ™ as noo [n = 0, ±1, ±2, ...]. (8) 

The integer n (the winding number) is related to the Jacobian of an S3 — > S3 
map and is given by [H] 

n= Vh?j €i ^ A n A n A n- ( 9 ) 

This expression is closely related to the Bardeen current K^. Indeed, in the 
A° = gauge only K° 7^ and one finds for pure gauge fields: 

K° = -^ k e abc A' a AiA k c = ^ige ljk Tr A i A j A k . (10) 

The true vacuum is a superposition of these, so-called, n-vacua and is 
called the 0-vacuum: 

\6 >= S n e- in9 \n > . (11) 

It is easy to see that in the vacuum to vacuum transition amplitude there are 
transitions with j da^K^ 7^ 0. Indeed, 

n\ t=+00 - = ^ J da^\$Z±™. (12) 

Using Eq (12), in detail one can write for the vacuum to vacuum transition 
amplitude 

+ < e\6 >_ = S m n e m9 e- m9 + < m\n >_ = S v e iv8 S n + < n + v\n >_ . 

(13) 

It is easy to see that the difference in winding numbers v is given by 

" = 32^ / ^"IKS - jj^ / A^rV (14) 

Using the usual path integral representation for the vacuum to vacuum am- 
plitude _|_ < 6\6 >_, one sees that 
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+ < 0\6 >_= K J 5Ae lS ^5[v - gf^ / dtxFFFap,], (15) 

where 

S cff L4] = S QC d[A] + 9^_ J d A xF^F ailv . (16) 

The resolution of U(1)a problem, by recognizing the complicated nature of 
the QCD's vacuum, effectively adds and extra term to the QCD Lagrangian: 

This term violates Parity and Time reversal invariance, but conserves Charge 
conjugation invariance, so it violates CP. The existing strong bound on the 
neutron electric dipole moment [7j \d n \ < 3 x 10 ~ 26 ecm requires the angle 9 
to be very small [d n ~ e9m q /Mf i implies [HJE] 6 < 10 -9 ]. Why should this 
be so is known as the strong CP problem. 

This problem is actually worse if one considers the effect of chiral trans- 
formations on the 9- vacuum. Chiral transformations, because of the anomaly, 
actually can change the 9- vacuum [TU| : 

e iaQs \6>=\e + a> . (18) 

If, besides QCD, one includes the weak interactions, the quark mass matrix 
is in general complex: 

^Mass = qiRM i:j q jL + h.c. (19) 

To go to a physical basis one must diagonalize this mass matrix and when one 
does so, in general, one performs a chiral transformation which changes 9 by 
Arg DetM. So, in the total theory, the coefficient of the FF term is 

9 = 9 + Arg detM (20) 

The strong CP problem is why is this 9 angle, coming from the strong and 
weak interactions, so small? 



2 Approaches to the Strong CP Problem 

There are three possible "solutions" to the strong CP Problem: 

i. Unconventional dynamics 

ii. Spontaneously broken CP 

iii. An additional chiral symmetry 

However, in my opinion, only the third of these is a viable solution. Of course, 
it might be possible that, as a result of some anthropic reasons 9 just turns out 
to be of O(10 -10 ) but I doubt it, as a Universe where CP is violated strongly 
seems as viable as one where it is not. 
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Appealing to unconventional dynamics to solve the strong CP problem 
is also not very believable. Attempts have been made to suggest that the 
boundary conditions which give rise to the 9- vacuum are an artifact |12) . but 
then what is the solution to the U(1)a problem? Other approaches try to use 
the periodicity of vacuum energy E{9) ~ cosO to deduce that 9 vanishes |T5|. 
but fail to motivate why one should minimize the vacuum energy. 

The second possibility, of spontaneously broken CP, is more interesting. 
Obviously, if CP is a symmetry of nature which is spontaneously broken, then 
one can set 9 = at the Lagrangian level. However, if CP is spontaneously 
broken 9 gets induced back at the loop-level. To get 9 < 10 -9 one needs, in 
general, to insure that 9 vanishes also at the 1-loop level. Although models 
exist where this is accomplished Jl] , theories with spontaneously broken CP 
need complex Higgs VEVs, leading to difficulties with FCNC and domain 
walls an d one must introduce recondite physics ^3 E| to avoid these 
problems. In my view, however, the biggest drawback for this "solution" to 
the strong CP problem is that experimental data is in excellent agreement 
with the CKM Model- a model where CP is explicitly, not spontaneously 
broken. 

Introducing an additional chiral symmetry is a very natural solution for 
the strong CP problem since this chiral symmetry, effectively, rotates the 9 
-vacua away. Two suggestions have been made for this chiral symmetry: 

i. The u-quark has no mass, m u — ^B] 

ii. The Standard Model has an additional global U(l) chiral symmetry ^2] 
The first possibility is disfavored by a standard current algebra analysis [20] 
which shows that all data is consistent with a finite mass for the u-quark. 
Furthermore, it is difficult to understand what would be the origin of this 
chiral symmetry, which asks effectively that Arg detM = 0. 

3 U(1)pq and Axions 

Introducing a global chiral U(l) symmetry |19] - which has become known 
as a U(1)pq symmetry- provides perhaps the most cogent solution to the 
strong CP problem. This symmetry is necessarily spontaneously broken, and 
its introduction into the theory effectively replaces the static CP-violating 
angle 8 with a dynamical CP- conserving field- the axion. The axion is the 
Nambu- Goldstone boson of the broken U(1)pq symmetry. j2J As a result, 
under a U(1)pq transformation, the axion field a(x) translates 

a(x) -» a(x) + af a , (21) 

where f a is the order parameter associated with the breaking of U(1)pq. 

Formally, to make the Lagrangian of the Standard Model U{\) pq invariant 
this Lagrangian must be augmented by axion interactions: 
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itotai = L SM +6^F^F a ^-\d^a+L ini [d^a/f a -,^H^i^K V Fa^- 
OZ7T z j a 6Zn 

(22) 

The last term above is needed to ensure that the U(l) pq current indeed has 
a chiral anomaly: 

dfiJpQ = ^ 3271-2 Pafj,u- (23) 

This term also represent an effective potential for the axion field, and its 
minimum occurs at < a >= — 

Since at the minimun the 0-term is cancelled out, this provides a dynamical 
solution to the strong CP problem. JHj 

It is easy to understand the physics of the Peccei Quinn solution to the 
strong CP problem. If one neglects the effects of QCD then the extra U(\)pq 
symmetry introduced allows all values for < a > to exist: < < a > < 2ir. 
However, including the effects of the QCD anomaly serves to generates a 
potential for the axion field which is periodic in the effective vacuum angle 

e + i^-. 

V e& ~cos[e + i^±}. (25) 

J a 

Minimizing this potential with respect to < a > gives the PQ solution: 

< a >= -I±0, (26) 

Obviously, the Lagrangian (22) written in terms of a p h ys = a— < a > no 
longer has a CP violating 0-term. 

Expanding V e s at the minimum gives the axion a mass 

2 d 2 V cS ±_^_d_ » v p . , 

The calculation of the axion mass was first done explicitly by current algebra 
techniques by Bardeen and Tye. [22] Here I will give an effective Lagrangian 
derivation of m a , as this technique also gives readily the couplings of 
axions to matter. 



4 Axion Dynamics 

In the original Peccei Quinn model, the U{\)pq symmetry breakdown coin- 
cided with that of electroweak breaking /„ = vp, with vf ~250 GeV. However, 
this is not necessary. If f a >> vp then the axion is very light, very weakly 
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coupled and very long lived. Models where this occurs have become known as 
invisible axion models. 

It is useful to derive first the properties of weak-scale axions and then 
generalize the discussion to invisible axion models. To make the Standard 
Model (SM) invariant under a U(1)pq transformation one must introduce 2 
Higgs fields to absorb independent chiral transformations of the u- and d- 
quarks (and leptons). The relevant Yukawa interactions involving these Higgs 
fields in the SM are: 

-^Yukawa — 



Defining x — v^jvi and vp = y/vf + v%, the axion is the common phase field 
in and #2 which is orthogonal to the weak hypercharge: 

$1 = ^e iax/vF [l] ; <Z> 2 = ^=e la/xVF [\ (29) 
v2 " v 2 1 

From the above, it is clear that ^Yukawa is invariant under the U(1)pq trans- 
formation 

a^a + avp; u Rj - e~ lax u R] ; d Rj -» e~ ia ' x d Rj ; £ Rj - e~ la ' x l Ry (30) 
Let us focus on the quark pieces. The symmetry current for U(1)pq, 

Jp Q = -vpd^a + xSiUm^um + -HidiR^diR (31) 

identifies the anomaly coefficient £ in Eq. (23) as: 

N 

£=-{x + l/x)=N g {x + l/x). (32) 

To compute the axion mass, it is useful to separate the effects of the axion 
interactions with the light quarks from the rest. These interactions can be 
deduced from the underlying theory by constructing an appropriate effective 
chiral Lagrangian. The effects of the heavy quarks essentially can then be 
accounted through their contribution to the chiral anomaly of Jpq 

For the two light quarks, as usual |5U], one introduces a 2 x 2 matrix of 
Nambu-Goldstone fields 

£ = ex P [i^±^], (33) 

where f„ is the pion decay constant. Then the meson sector of the light quark 
theory, neglecting the effect of the Yukawa interactions, is embodied in the 
U(2)v x U(2)a invariant effective Lagrangian 

Chiral = -ffTrd^E^. (34) 
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To ichirai one must add U(2)v x U(2)a breaking terms which mimic the 
U{\)pq invariant Yukawa interactions of the u- and d-quarks. This is accom- 
plished by the Lagrangian 

L masa = \{UmlfTr[EAM + (EAM)% (35) 

where 

-iax/v F n Hhi f) 

A=f n y, XVF ] ; M m „. ]. (36) 

m u +m d 

Note that the invariance of L mass under Z7(l) pq requires that £ transform as 

iax n 

S ^Z[ e wJ- (37) 

L mass , however, only gives part of the physics associated with the sym- 
metry breakdown of U{2)a- In fact, the quadratic terms in L mass involving 
neutral fields 

r(2) ( m Tr) 2 f m u r o , x fir i2 , m d r o fir i2i / O o\ 

mass 2 L TO u + TO d L w F J m u + m d u xv p 1 ] y ' 

give for the ratio 

^ = 24 K 1.6. (39) 

which contradicts experiment. Indeed, if L mass was all that there was, we 
would have recreated in the effective Lagrangian language the U(1)a problem! 
Furthermore, with only this term the axion is still massless. 

The resolution of the U(1)a problem in the effective Lagrangian theory 
is achieved by adding a further mass term which takes into account of the 
anomaly in both U(1)a and U(X)pq. This mass term gives the r\ the right 
mass and produces a mass for the axion. It is easy to convince oneself that 
such a term has the form [231 



K) 2 r | fir (N g -l)(x + l/x) 12 

^anomaly = 7T~ IV H ~ n G > ( 40 

' Vp 



where (m°) 2 ~ m 2 >> m 2 . The coefficient in front of the axion field in 
-^anomaly reflects the relative strength of the couplings of the axion and the r\ 
to FF as the result of the anomalies in U(1)pq and U(1)a- Naively, the ratio 
of these couplings is just f 7r /2vF^. However, the reason that N g — 1 appears 
in the above, rather than N g , is that L mass already includes the light quark 
interactions of axions, so only the contribution of heavy quarks to the PQ 
anomaly should be taken into account in L anoma i y . 

Diagonalization of the quadratic terms in L mass and £ anoma i y gives both 
the axion mass and the parameters for axion-pion and axion-eta mixing for 
the PQ model. It is convenient to define 
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25 keV. (41) 



Then one finds: 



where 



m a = X m m a ; £, av — A 3 — ;£ ai) = A — , (42) 



A m = N g (x + l/x) (43) 
A 3 = h(x- l/x) N g (x + l/x)?±^±] (44) 

\ = ±(l-N g )(x + l/x). (45) 

In addition to the three parameters above, all axion models are character- 
ized also by how the axion couples to two photons. Writing the interaction 
Lagrangian describing this coupling as 



L = —K flphys ' F^F (46) 

47T )a 



one needs to find the coupling K ail iai the PQ model. This coupling follows 
from the electromagnetic anomaly of the PQ current 

«q = ^F^. (47) 

Here £ 7 gets contributions from both quarks and leptons, and one finds: 



£ 7 = JV fl {[3(2/3) a ]a; + [3(-l/3) 2 + (-1) 2 ]!/*} = -N g (x + l/x). (48) 



As before, in computing K ail one must separate out the light quark contribu- 
tion of the axion in the anomaly ( so that £° ff - = ^N g (x + l/x) — \x— |A) and 
add back the axion to two-photon contribution which arises from the coupling 
of the ir° and r\ to two photons, via axion-pion and axion-eta mixing: A3 + IA0. 
This then gives 

K ail = N g {x + l/x)[ (49) 

m u + m d ' 



5 Invisible Axion Models 

The original PQ model jlH) . where f a = vf, was long ago ruled out by exper- 
iment. For example, one can estimate the branching ratio |23| 

BR(K + -> vr+ + a) ~ 3 x 1(T 5 A§ = 3 x !Q- 5 (x + l/x) 2 (50) 



10 R.D. Peccei 



which is well above the KEK bound £3] BR(K+ — > tt+ + nothing) < 3.8 x 
10 -8 . However, invisible axion models, where f a » vp, are still viable. 

Invisible axion models introduce scalar fields which carry PQ charge but 
are ST/ (2) x U(l) singlets. As a result, it is possible that the VEVs of these 
fields have a scale much larger than the one set by the weak interactions. 
Basically, two types of models have been proposed. The first of these models, 
due to Kim j2S] and Shifman, Vainshtein and Zakharov [2HJ (the, so-called, 
KSVZ Model), introduces a scalar field a with f a =< a > >> vf and a 
superheavy quark Q with Mq ~ f a as the only fields carrying PQ charge. The 
second of these models, due to Dine, Fischler and Srednicki and Zhitnisky 
28 (the, so-called, DFSZ Model), adds to the original PQ model a scalar field 
<fi which carries PQ charge and f a =< <j) > >> Dp. 

It is straightforward to repeat the calculations we just did above for these 
models, to get the axion mass and couplings. I will do this explicitly here 
for the KSVZ model, because it is simple and illustrates well what we just 
did, and will just quote the results for the DFSZ model. By construction, the 
KSVZ axion does not interact with leptons and it only interacts with light 
quarks as the result of the strong and electromagnetic anomalies 

where eQ is the em charge of the superheavy quark Q. 

Since in the KSVZ model the ordinary Higgs do not carry PQ charge, 
the only interactions of the axion with the light quark sector come from the 
effective anomaly mass term, which here is given by 

^anomaly = ~ ~ ^ [V + Trr®] 2 (52) 
* *Ja 

To the above one must add the standard quadratic term coming from the light 
quarks 

4 2 L = -^i^r(-° + vf + -^-{n - o 9 ]. (53) 

2 m u + m d m u + m d 

(2) 

Diagonalizing L an0 maiy and Lm« gives 

Vf _ , vf _ t \ /t c \ U 
m a = —m a = \ m —m a ; £ a7r = A 3 — ; t, arl = A — , (54) 

J a J a J a J a 

where for the KSVZ model these parameters are: 

Am = 1 ! A 3 = -— ■ r ; A = --. (55) 

2(m u + m d ) 2 



Note that, numerically, in the KSVZ model the axion mass is given by the 
formula: 
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m „ = ^e.3[^]eV. (56) 

J a J a 

The calculation of K ail in this model is equally straightforward. To the 
contribution of the superheavy quark in the electromagnetic anomaly [3eg] , 
one must add that coming from the mixing of the axion with the ir° and the 
V [ A3 + §Aq ]. This gives, finally, 



2 4m d + m u 
3(m u + md) 



Ka~fy — 3eg — r)/ — j ^ . (57) 



I will not go through the analogous calculation for the DFSZ model, but 
just quote the results. For this model, it proves convenient to define 

_ 2w| 2v\ 
X\ — — =- X2 = (58) 
vf v F 



where vf = \Jv\ + v\ and v\ and vi are the two Higgs VEVs. Furthermore, 
if one rescales f a — > the axion mass in the DFSZ model is given by the 
same equation as for the KSVZ model, Eq. (56), corresponding to A m = 1. 
With this rescaling understood, one then finds: 

1 Xi-X 2 m d ~m u 1 — iV s . . 

A3 = O ^ 1 5 A = OAT ' 59 

2 2N g md + m u 2J\ g 



and 

4 _ 4m rf 
3 3(m d + m u )' 

Although the KSVZ and DFSZ axions are very light, very weakly coupled 
and very long-lived, they are not totally invisible. I will not describe in detail 
how these "invisible" axions affect astrophysics and cosmology, as others in 
this meeting have done so. However, in an attempt to make this brief precis on 
the strong CP problem and axions self contained, I will add a few comments 
on the bounds one derives for invisible axions. 

Astrophysics gives bounds on the axion mass since axion emission, through 
Compton production [ye — > ae] and the Primakoff process, causes energy 
loss in stars. The energy loss is inversely proportional to and hence 
proportional to m\. Thus axions must be light enough, so as not to affect 
stellar evolution. This is not the only way in which one obtains an upper 
bound on m a . For instance, another bound on m a come from SN1987a, since 
axion emission through the process N N — ► NNa in the core collapse affects 
the neutrino spectrum. |3(J) Typical bounds |31j obtained from astrophysics 
require axions to be lighter than 

m a < 1 - l(T 3 eV. (61) 

Remarkably, cosmology gives a lower bound on the axion mass (and an 
upper bound on f a ). |32l 1331 The physics connected with this bound is 
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simple to understand. When the Universe goes through the PQ phase tran- 
sition at T ~ /„ >> Aqcd the QCD anomaly is ineffective. As a result, 
early in the Universe < a p h ys > is arbitrary. Eventually, when the Universe 
cools down to temperatures of order T ~ Aqc d the axion obtains a mass and 
< a P hys >^ 0. The PQ mechanism through which this happens, however, is 
not an instantaneous process and < a p h ys > oscillates to its final value. These 
coherent p a = axion oscillations contribute to the Universe's energy density 
and axions act as cold dark matter. The energy density contained in axion 
oscillations is proportional to f a and thus bounds on the energy density of 
cold dark matter in the Universe provide an upper bound on f a (and a lower 
bound on m a ). 

I quote below the result of a recent calculation by Fox, Pierce and Thomas 
13*5] of the contribution to the Universe's energy density due to axions: 

Qah2 = °- 5[ iMev ]i[e * + m2]t (62) 

Here £ is the coefficient of the PQ anomaly and, with the way we defined f a , 
£ = 1 for both the KSVZ and DFSZ models. The angle 6{ is the misalignment 
value for < a p h ys > /fa and ag is its mean squared fluctuation. Finally, 7 is a 
possible dilution factor for the energy density produced by axion oscillations. 

One can use the recent WMAP bound on cold dark matter to bound 

n a h 2 < 0.12 (63) 
Assuming no dilution (7 = 1) and using for 9i an average misalignment angle 

rj 2 

9f = 2jp and neglecting its fluctuations, the WMAP data gives the following 
cosmological bound for the PQ scale: 

f a < 3 x 10 n GeV or m a > 2.1 x !CT 5 eV. (64) 



6 Concluding Remarks 

My conclusions are rather simple and telegraphic. In my view, after more than 
25 years, the preferred solution to the strong CP problem remains still the 
idea that the Standard Model has an additional U(1)pq symmetry. Such a 
solution, necessarily, predicts the existence of a concomitant axion. Although 
Fermi scale axions have been ruled out, invisible axions models are still viable 
and axion oscillations toward its minimun could account for the dark matter 
in the Universe. This is an exciting possibility. However, no totally compelling 
invisible axion models exist, and there are no strong arguments to believe that 
f a takes precisely the value needed for axions to be the dark matter in the 
Universe. Nevertheless, it is encouraging that experimentalists are actively 
searching for signals of invisible axions. 
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